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' The invariant mass spectra of the n'^n and tt^tt" pairs in the 

process e'^e^ tv^tv'tv'^ were studied in the SND experiment at the 
VEPP-2M colhder in the energy region^ ~ 1020 MeV. These studies 
were based on about 0.5 x 10® experimental events. The spectra were 
^ , analyzed in the framework of the vector meson dominance model. It 

was found that the experimental data can be described with e e — + 
>J I pTT — > 7r^n~-K^ transition only. Upper limit on the branching ratio 

■ of the <?!>(1020) Tr^vr^Tr" decay through intermediate states different 

from pvr was obtained at the 90 % confidence level: B{(f> — > 7r"'"7r~7r'') < 
6- 10"''. The p-meson mass and width which follow from the spectra 
analysis are rUp = 775.0 ± 1.3 MeV, Fp = 150.4 ± 3.0 MeV. Neutral 
and charged p-mesons mass difference was found to equal mp± — 771^0 = 
^ , —1.3 ± 2.3 MeV. In the tt^tt" mass spectrum the p — uj interference 

was seen at two standard deviations level. 

Oh! 

^ ■ I. INTRODUCTION 
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In the framework of the vector meson dominance model (VDM) the cross section 
of the process e^e" — *■ tt+tt^tt'^ in the energy region ~ 1 GeV is determined 
' by the amplitudes of vector mesons Y {V — uj,(j),uj' ,uj") transitions into the final 

state: V tt+tt^tt^. From experimental data it is known that the pir intermediate 
state dominates in these transitions. In the energy region ^/s ~ the main 
contribution to the cross section comes from the (^(1020)-meson decay (p ^ pn ^ 
7r^7r~7r°. The total cross section of the process e^e~ 'W'ff~'''° vicinity 
of the (/)(1020) resonance was studied in several experimentsErEl, including the SND 
(Spherical Neutral Detector) experiment at VEPP-2M coUideiQ. Studies of tiiis 
process dynamics (the dipion mass distribution) were reported only in two worka3i3. 

The transition V 7r+7r~7r° besides pn intermediate state (Figj^a) can be per- 
formed by p'^"^TT intermediate state (Figj^c). Even if the branching ratio of the 

V p'^"^n decay is rather small, the interference of this amplitude with the 

V ^ pn amplitude can give the noticeable contribution to the total cross sec- 
tion. The transition e+e" —^ Tr+Tr^Tr*^ is also possible through p — uj interference: 
e+e" y W7r° ^ p°t:^ {V = p,p',p") (Fig.|l|b). This effect was observed 
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in the SND experiment in the energy region -/s = 1200-1400 MeVtl. In the pro- 
cess e~^e'^—^ pTT the interaction of the p- meson with pion in the final state is also 
expectedN (Fig||d). 

Besides of the — > Stt transition dynamics studies, the process e+e^ tt+tt^tt^ 
in the 0(1020) resonance energy region can be used for the p-meson parameters 
measurement. In the reaction e^e~ pi: the neutral and charged p-metpiis are 
produced, so it can be used for the and mass difference determinatior£iL-|The 
p mass values obtained by using different reactions contradict to each otherEHI. So 
it is worthwhile to compare the resonance parameters {nip and Fp), obtained from 
the e+e^ — )■ pi: reaction with the results of other experiments. 

Investigation of the tttt mass distribution in the process e'^e~ —> 7r"*'7r~7r'^ provides 
important information about vector mesons and their interference. Here we present 
the results of the dipion mass spectra analysis in SND experiment. 

II. DATA ANALYSIS 
A. Experiment 

The SND detectorEl operated since 1995 up to 2000 at the VEPP-2M0 colfider 
in the energy range-^/s from 360 to 1400 MeV. The detector contains several subsys- 
tems. The tracking system includes two cylindrical drift chambers. The three-layer 
spherical electromagnetic calorimeter is based on Nal(Tl) crystals. The muon/veto 
system consists of plastic scintillation counters and two layers of streamer tubes. 
The calorimeter energy and angular resolution depend on the photon energy as 
oe/E{7o) = 4.2%/^£;(GeV) and a^^.g = QM° /^E{GeV) ® 0.63°. The tracking 
system angular resolution is about 0.5° and 2° for azimuthal and polar angles re- 
spectively. The energy loss resolution dE/dx in the drift chamber is about 30% 

- good enough to provide charged kaon identification in the (/)-meson production 
region. „ 

In 1998 SND had collected data in the energy region from 984 to 1060 MeVtl with 
integrated luminosity about 8.5 pb"""^. For studies of the dipion mass spectra in the 
process e~^e~ 7r+7r~7r° the data sample collected in the energy region-^s = 1016 

- 1022 MeV was used. The total integrated luminosity accumulated in this region 
is about 4.3 pb^^. 

B. Selection of e^e~ — > tt'^tv^tt'^ events 

The data analysis and selection criteria used in this work are similar|_to those 
described in Ref.El. During the experimental runs, the first-level triggeiO selects 
events with energy deposition in the calorimeter more then 200 MeV and containing 
two or more charged particles,|-puring processing of the experimental data the event 
reconstruction is performedEji]. For further analysis events containing two or more 
photons and two charged particles with \z\ < 10 cm and r < 1 cm were selected. 
Here z is the coordinate of the charged particle production point along the beam 
axis (the longitudinal size of the interaction region is cr^ ^ 2.5 cm); r is the distance 
between the charged particle track and the beam axis in the r — (p plane. Extra 
photons in e+e~ — > tt'^tt^tt'^ events can appear because of the beam background 
overlap or nuclear interactions of the charged pions in the calorimeter. Under these 
selection conditions, the background sources are e"^e~ — > K^K~ , KsKl, rj^, a-iTr*^, 
e"'"e~77 processes and the beam background. 



2 



To suppress the beam background, the following cuts on the angle between the 
two charged particles tracks ip and energy deposition of the neutral particles Eneu 
were applied: ip > 40°, Eneu > 0.1 ■^/s. 

To reject the background from the e"*"e^ K^K^ process the following cuts 
were imposed: (dE/dx) < 5 • {dE / dx)min for each charged particle, (dE/dx) < 
3- {dE/dx)min at least for one of them, and Acf) < 10°. Here |A(/)| is an acollinearity 
angle in the azimuthal plane and {dE/dx)„iin is an average energy loss of a minimum 
ionizing particle. To suppress the e^e^ — > e~^e~jj events the energy deposition of 
the charged particles Echa was required to be small enough: Echa < 0.5 -^/s. 

The selected events were reconstructed under assumption that all registered par- 
ticles have the common vertex in the e~^e~ interaction region. The average posi- 
tion of the beams interaction point and its rms spread were measured by using 
e'^e~ e^e~ events during data collection. Then a kinematic fit was performed 
under the following constraints: the charged particles are considered to be pions, the 
system has zero total momentum, the total energy is-\/s, and the photons originate 
from the 7r° 77 decays. The value of the likelihood function is calculated 
during the fit. In events with more than two photons, extra photons are considered 
as spurious ones and rejected. To do this, all possible subsets of two photons were 
inspected and the one, corresponding to the maximum likelihood was selected. The 
two-photon invariant mass and xi^ distributions are shown in Fig.|^ and Fig.^. Af- 
ter the kinematic fit the following additional cuts were applied: 36° < 6^ < 144°, 
= 2, and Xstt < 20. Here 9^ is polar angle of one of the photons selected by 
the reconstruction program as originated from the Tr'^-decay, is the number of 
detected photons. Under these criteria about 0.5 x 10^ events were selected. In 
Fig.^ and Fig.|| the angular distributions of pions for the selected events are shown. 
While Fig.^ and Fig.^^ demonstrates the photon energy distributions for the same 
events. _ 

The difference between described selection criteria and those used in Refill is that 
the cuts on the pion and photon angles were relaxed here. This leads to the increase 
of the selected events number by a factor of about 1.5. The (j) meson parameters 
obtained under this new selection criteria agree with results reported in the previous 
work. For example, the e+e^ tt+tt^tt^ cross section changed by about 0-5%, what 
is negligible in comparison with the 5% systematic error declared in Ref.LJ. 

The number of background events was estimated by using simulation in the fol- 
lowing way: 

Nbkgis) = C^ffi(s)e^(s)-^i(s), (1) 

i 

where i is a process number, aui{s) is the cross section of the background process 
taking into account the radiative corrections, IL{s) is the integrated luminosity, 
ei(s) is the detection probability for the background process. The accuracy of the 
background events number determination is estimated to be about 15 - 20 %. The 
numbers of e+e~ — > tt^tt^tt'' events (after background subtraction) and background 
event numbers are shown in Table H. The e^e~ — > KsK^ events are the main 
background source and their contribution is ^ 1% of all selected events. 

C. Construction of the tt^tt^ and tt^tt" mass spectra 

For the selected events the tt+tt" and tt^tt" mass spectra were formed and ar- 
ranged in histograms with dipion mass range from 280 to 880 MeV and bin width 
of 20 MeV. All experimental energy points were used to produce a single invariant 
mass distribution. In case of the tt^tt" distribution the charged pion was selected at 
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random from the two possibilities. The invariant mass values were calculated after 
the kinematic reconstruction. 

The expected background was subtracted bin by bin while forming the desired his- 
tograms. Simulated distributions were used for subtraction of the e'^e^ K^K^, 
rij, lott'^ and e~'"e~77 backgrounds. The e~^e~ —> KsKl main background was 
studied using experimental events. 

As a rule, the background from the e+e^ KsKl originates when Ks 
Tr~^TT~ decay inside the collider vacuum chamber is accompanied by "photons" in 
calorimeter produced due to Kl meson nuclear interactions or its decay in flight 
inside the calorimeter. To select KsKl events for related background studies the 
following additional cuts were applied r > 0.2 cm for both charged particles and 
-0 > 130°. Note that Ks decay length is about 0.5 cm for our energies and the 
angle ij} between the pions from the Ks decay is about 150°. 

The experimental and simulated "tt^tt'^" distributions are consistent to each other 
(Fig|8|), so the simulated distribution was used for the KsKl background subtrac- 
tion in the tt^tt^ invariant mass spectrum. The agreement is not so good for the 
"tt+tt"" mass spectrum: the simulated histogram is shifted from the experimental 
one by about 3% in the peak region (Fig|9|). This disagreement is caused due to 
inaccuracies in simulation of the K^ nuclear interactions. The experimental distri- 
bution (Fig.^ was used for KsKl background subtraction for events in which both 
charged particles had r > 0.2 cm. Unfortunately, this procedure is not justified if 
at least one charged particle in the event has r < 0.2 cm. Indeed, "tt+tt"" mass 
distribution for the KsKl background depends crucially on the shape of the angle 

distribution, which in its turn depends on the value of r, because during event 
reconstruction it was assumed that all charged particle tracks have a common ver- 
tex in the beams interaction region. For the Ks tt+tt" vertex this assumption 
is not true. Therefore, "tt+tt"" distribution obtained for KsKl background un- 
der condition r > 0.2 cm cannot be used for subtraction of all KsKl background. 
Instead for events with r < 0.2 cm for one of the charged particles the simulated 
distribution, corrected for the above mentioned 3% discrepancy, was used for KsK^ 
background subtraction. The dipion invariant mass spectra and background con- 
tributions obtained in such a way are shown in Fig.|lO|, O. 



III. THEORETICAL FRAMEWORK 



In the VDM framework the cross section of the e+t 



TT+vr tt'^ process is 
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where and p_ are the 7r+ and tt momenta, too and to+ are tt+tt and tt+tt*^ 
invariant masses. Formfactor F of the 7* —^ Tr+vr^vr*' transition has the form 
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where m_ is 7r~7r'^ invariant mass, m-^o and ttt-tt are the neutral and charged pion 
masses. The — s- tt^tt" and — > tt^tt" transition couphng constants could be 
determined in the following way: 

2 67rTO^^,rpO 2 eTrm^^Tpi 



If these constants are equal, then the and p^ meson widths are related as follows: 



"^p± 9o(,"ipoj'^ 

Factor Z{ni) = 1 — isi'^{xu,s) takes into account the interaction of the p and 
TT mesons in the final stateEj. <i>(m, s) as a function of the pion pair invariant 
mass m is shown in Fig.^ for the e^e" center-of-mass energy = m^. The free 
parameter in our fit si was introduced to account for possible deviations from the 
Ref.li3 jprediction si = 1. The fourth item in (||) takes into account the p — to 
mixingij. Polarization operator of these raiding Ifpi^ satisfies lm{Ilp^) ^ Re(np^), 



where Re(np^^) = 2m^6 with S = 2.3 MeVEJO. So we have assumed Im(np^^) = 
in subsequent analysis. Amplitudes of the 7* —^ pir and 7* —^ tjvr" transitions have 
the form 

V = LO,<p,LO',<^" ^ ^ v = p,p',p" ^ ^ 

Using SJjll) measurements of the e+e^ tt+tt^tt^ and e+e^ tt'^tt'^'^ cross 
sectionsQ'EII, we can express the ratio of these amplitudes in the 0-meson region 
as 

A^^{s)/Ap^{s) = K^T,/p^e'"^'-''^'"',K^^/pT, ~ 0.3,0^^^/p^ ~ -110°,v^ = m^. 

The amplitude aa^r takes into account possible additional intermediate states in the 
7* 7r~'"7r~7r'' transition. For example, 7* — s- p'{p")Ti — > tt^tt~tt^ transition will 
lead to the contribution 

Ap'Trjs) I gp'TTir ffp'TTTT ffp'TTTT \ 

" Ap^{s) \Dp,{mo) Dp,{m+) ^ Dp,{m_) j ' 

where 

The ratio ^p'7rMp;r is expected to be real for energies -^/s ~ mjll where 

^p'tt ^ 94>p''^ 

ApTT Q<pp7^ 

So the imaginary part of (^) is negligible due to the large p' mass (mp/ ~ 1400 
MeV). Therefore the a3,r amplitude was assumed to be real in our analysis. 
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IV. APPROXIMATION OF THE tttt MASS SPECTRA 



A. Theoretical distributions 



The experimental dipion mass spectra (Tables O and |T|) were fitted with theo- 
retical distributions. Using the e+e^ 7r+7r~7r*'cross section (^) the theoretical 
spectra were calculated in the following way: 



jTij+A m™''^(mo) 
mj-A m™"'(mo) 



^r(-'*) = 7^- / ""'odrno / m+dm+\p+ x p_\^ ■ \F\^ , (8) 



^fls) = ^^- / m±dm± I modmo\p+xp_\' -IFI', (9) 



-A 



where j is the bin number, A = 10 MeV - half of the bin width, nij - the central 
value of the invariant mass in the jth bin, C's{s) - normalizing coefhcient. These 
spectra were corrected taking into account the detection cfBcicncy ej (Fig.|l3| and |l^ 
for the jth bin and a probability for the event belonging to the jth bin to migrate 
to the ith bin due to finite detector resolution 

g["\s) = o^sfis)^^ . (1 + Sr\s)) . /3f (10) 



Here 6i{s) is a radiative correction, /S^-"^ and /3|^^ are corrections due to inaccuracies 
in simulation of the nuclear interactions of the charged pions. The values of a^j, ej 
and Si{s) were obtained from simulation. 

These distributions were formed for all energy points ^/s (Table |) and then 
summed as follows: 

Here Cc{s) = '^iGi{s) is a normalization factor, w{s) = N{s)/N is the weight 
factor, N{s) and N being the numbers of the e+e^ 7r+7r~7r° events in each 
energy point and the total events number respectively (Table |). 



B. Detection efficiency 

To obtain the Uij, ej and 5i{s), the Monte-Carlo sample of about 1.4x 10^ e+e^ — > 
pir —^ tt+tt^tt'' events were processed in the same way as experimental data. The 
statistical errors of a,jj coefficients were included in the error of P, (Tables.0 and 
[H). The average detection efficiency obtained by simulation is about 0.37 and 
varies from 0.05 to 0.4 depending on dipion mass value (Fig. |l|and 0). 
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Inaccuracies, in the Xstt ^^^^ dE/dx simulations lead to an error in average detec- 
tion efRciencyQ, but do not modify its dependence on the tttt invariant mass, except 
the regions of m^rTr < 315 MeV and m^^ > 845 MeV. The cut iV^ = 2 is also a 
source of some error in the detection efficiency determination, due to inaccuracies in 
the extra photons simulation. As was mentioned above, extra photons can appear 
because of the beam background overlap or charged pion nuclear interactions. In 
fact, if the extra photons originate from the beam background overlap, — 2 
criteria does not modify the detection efficiency dependence on a dipion mass, be- 
cause the probability of beam background overlap does not depend on dipion mass. 
On the contrary, the probability of extra photons appearance because of nuclear 
interactions of the charged pious does depend on pion energy. Therefore any inac- 
curacy in simulation of these nuclear interactions will transform in inaccuracies in 
the simulation of extra photons and can lead to the detection efficiency error which 
will depend on the dipion mass. To correct this error we used the correction factor 
obtained in the following way: 

where Ki and Ni are the event numbers in the ith bin of experimental and simulated 
mass spectra respectively under condition N-^ > 2. ki and rii are the event numbers 
when the cut = 2 was imposed. The dependence of these correction factors on 
the dipion mass (Fig.p^ can be approximated by linear functions: 

/jf) = • (1 + ao • m,), (3^1^^ = C(±) • (1 + a± • m,) 

with slopes oq = —0.13 • 10^"^ and a± = 0.083 • 10"'^ for tt+tt^ and tt^tt^ pairs 
respectively. As was mentioned above, the detection efficiency in the regions TOttu- < 
315 MeV and 

ttZtttt ^ 845 MeV have a large error. Uncertenity of Pi in these regions 
was increased by addition of the 100% spread. 



C. Fitting of the experimental distributions 

The TT+TT^ and tt^tt^ mass spectra were fitted together. The function to be 
minimized was = Xq + X± i where 




Here H'-"^ and H'-^^ are the normalized tt+tt and tt^tt" mass distributions (his- 
tograms); erf ^ = Afff ^ ® AP^ and cr|=^^ = Ai?^^^^ APj;^^ include the uncer- 
tainties AHi and APi of the experimental and theoretical distributions. 

The fitting was performed with rUpO , mp± — WpO , TpO , Tp± as free parameters 
under the following assumptions: 

1. = 0, K^T,/p^ = 0, si = 0; 

2. was free parameter, K^^/p^,. = 0, si — 0; 

3. a^TT — 0, K^T^/pT^ and 4>^.^/pTr were free parameters, si = 0; 

4. 03^ = 0, K^^/p^ = 0, Si was free parameter. 
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The value for all variants approximately equals to 73 for 54 ~ 52 degrees of free- 
dom, besides Xo — 44, X± — 29. The difference Hi — Pi between the experimental 
and theoretic spectra are shown in Fig.|l^ and O. The difference for the tt+tt" 
distribution has a systematic spread of about 1%7 This spread was added to un- 
certainty of the theoretical P/^"* distribution and the fit was redone (Table.|^). For 
all variants of fitting the values of p meson widths TpO and Tp± are consistent with 
relation (||), and fits were performed again under assumption that gpOj^^ = gp±TrTi 
(Table.0) (r po was free parameter and r„± was calculated using expression (M)). 
The and p^ mass difference is consistent with zero, and the fits were redone 
again under assumption rripi = mpo (Table [v| ) . The experimental spectra fitted 
with theoretical distributions (variant 3 in Table |v^ ) are shown in Fig.^ and 

The parameters determined by fit are subject of various systematic errors origi- 
nated from model dependence, from inaccuracy in background subtraction and in 
detection efficiency determination. To estimate errors due to uncertainty in back- 
ground subtraction the experimental distributions were approximated with the sum 
of three histograms corresponding to the theoretical spectrum for the effect, KsKl 
background and expected summary background from sources other then KsK^. 
The background events fractions were taken as free parameters of the fit. The 
values of these fractions obtained from the fit differs by less than 20% from those 
estimated using equation (^. 

To study systematics related to inaccuracy in the extra photons simulation the 
corresponding correction factors were approximated by assuming the higher order 
polynomial, instead of linear function, to fit the dependence of these factors on the 
dipion mass shown in FigjT^. We also had tried the fit with linear approximation for 
the correction factors but with oq and a± as free parameters of the fit. The results 
of the fit ao = (-0.f3 ± 0.02) • lO'^ and a± = (0.12 ± 0.2) • lO'^ are consistent with 
those described above. 

The model dependence was estima ted by using results of the fits performed under 
different assumptions (Tables 0, 0, 0. 



V. DISCUSSION 



The SND results in comparison with other experimental data are shown in Table 



vn 



The fit results revealed that the experimental data can be described by using only 
/97r intermediate state. The value of additional, other than /ott, contribution a^T^ is 
consistent with zero (variant 2, in Table Vl): 



ag^ = (O.Of ± 0.23 ± 0.25) x 10"^ MeV^ 



-2 



Here the systematic error is related to uncertainties in background subtraction and 
detection efficiency determination. This result agrees with 03^ = (0.08 ± 0.45 ± 
0.37) X 10~^ MeV~^ reported in Re|U3. Using the a^-K value and results of the 
e^e'' TT+Tr^Tr" cross section studyQ the upper limit on the branching ratio of 
(f) Tr+TT^Tr*^ decay without the pir intermediate state can be obtained: 

B{4) vr+vr-vr") < 6 • lO"'^ (90% CL) 

This upper limit implies 

0.91 < — < 1.09 (90% CL) 

Here a^Ti is the total e+e^ —>■ tt^tt^tt^ cross section, CTpjr is e^e^ ~* pir ^ vt+tt^tt'^ 
cross section. The upper and lower limits above correspond to the constructive and 
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destructive interference between 03^ and pir related amplitudes. It may happen 
that e'^e" cj) p'lr amplitude is suppressed compared to e'^e^ ~f u^uj') — > p'lr 
amplitude. In this case one expects Im(a3,r) > R.e(a3T). So we had performed also 
the fit with Im(a3,r) as free parameter and supposing Re(a37r) = 0. The result of 
this fit also is consistent with zero: 

a37r = i- (-0.12 ±0.13 ±0.25) 



The fit with model taking into account p — lu mixing (variant 3 in Table VI ) gave 
the results 

K^^/prr = 0.20 ± 0.10 ± 0.05 = -125° ± 28°. 

The systematic error is related to uncertenity in background subtraction. The values 
obtained are consistent to the above given estimation -fCi^^/pyr — 0.3, 4>,^TT/pn — 
— 110°. The measured value of K^^jp^^ deviates from zero by about two standard 
deviations. 

Parameter si is related to the final state interaction of the p and tt mesons and 



was measured to be (variant 5 in Table VI) 



si = 0.3 ±0.3 ±0.3 

The main cause of the systematic error is again the background subtraction un- 
certenity. The result is consistent with zero, but also does not contradict to the 
prediction of Ref.EJ (si — 1 with about 10-20% accuracy). 

The measured neutral and charged p-mesons masses and their difference are (vari- 
ant 1 in Table 0): 

TOpO = 775.8 ±0.9 ±2.0 MeV, m^i = 774.5 ± 0.7 ± 1.5 MeV, 

mp± - rupa = -1.3 ± 1.1 ± 2.0 MeV. 

The systematic error of mp± is related to model dependence and uncertainties in 
detection efficiency determination. Background subtraction inaccuracy also con- 
tributes in the systematic errors of mpo and Jripi — mpo . The results obtained are 
consistent with ra^± = TOpO as well as with the world average for the mass difference 
—0.4 ± 0.8 MeVEJ. The measured mass-difference is also not in conflict with the 
prediction mp± — mpO = — 4.2± 1.2 MeVEd. If TOp± = WpO is assumed in the fit, the 



p-meson mass turns out to be (variant 1 from Table. VI) 



mp = 775.0 ± 0.6 ± 1.1 MeV 



This value is consi stent with the results of the e^e annihilation and r decay 



experiments (Table VII ). Note that, world average for these experiments is mp — 
776 ± 0.9 MeV. But the PDG valueB 769.3 ± 0.8 MeV, which takes into account all 
experiments in which the p-meson mass was measured, contradicts to the reported 
SND result. 



Our results for the neutral and charged p- meson widths are (variant 1 Table IV ) : 

Tpo = 151.1 ±2.6 ±3.0 MeV, Tpi = 149.9 ±2.3 ±2.0 MeV 

The systematic error includes model dependence, uncertainty in background sub- 
traction and in detection efficiency determination. Under assumption gpo-^^^ — .9p±7r7r 



(variant 1 Table VI) the p- meson widths were found to be 



Tpo = 149.8 ±2.2 ±2.0 MeV, Tpi = 150.9 ± 2.2 ± 2.0 MeV 



These results are consistent with other experipjental results (Table VII) as well as 
with the PDG world average 150.2 ± 0.8 MeVlB. 
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VI. CONCLUSIONS 



In the SND experiment at VEPP-2M collider in Novosibirsk the dipion mass 
spectra were studied in the e+e~ — > 7r+7r~7r'^ process at-^i = m^. The study is 
based on the data sample with about 0.5 x 10® experimental events. Spectra were 
analyzed within the VDM framework taking into account e~^e~ pn transition, 
p — u) mixing, final state interaction of the p and tt mesons, possible transition 
e~^e~ tt+tt^tt" through intermediate states different from pir (for example, via 
p'tt). Within the limits of our accuracy, the experimental data can be described 
with the e~^e~ pir 7r"'"7r~7r'^ transition only. The upper limit on the ^ Stt 
decay through intermediate states besides pn is B^cf) 7r+7r~7r°) < 6 • 10""^. In 
the 7r+7r~ invariant mass spectrum the p — ui interference was observed at two 
standard deviations level. The measured mass and width values for the neutral 
and charged p-mesons are nipo = 775.8 ± 0.9 ± 2.0 MeV, T^o = 151.1 ± 2.6 ± 3.0 
MeV, mp± = 774.5 ± 0.7 ± 1.5 MeV, T^i = 149.9 ± 2.3 ± 2.0 MeV. The difference 
between masses of the neutral and charged p-mesons was found to be mp± — rripO = 
— 1.3 ± 1.1 ± 2.0 MeV. Under assumption m^i = rUpO the following p-meson mass 
value was obtained rup = 775.0 ± 0.6 ± 1.1 MeV. Assuming the coupling constants 
equality 5pO^^ = (jpi^Trj the neutral and charged p-meson widths turn out to be 
r^o = 149.8 ± 2.2 ± 2.0 MeV and Tp± = 150.9 ± 2.2 ± 2.0 MeV. 
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(MeV) 




Nbkg 


NksKl 


Nk+k- 








~ 1017 


45913 


410 


331 


3 


42 


15 


19 


~ 1018 


101651 


1067 


901 


8 


104 


24 


29 


~ 1019 


137432 


1696 


1465 


13 


161 


26 


31 


~ 1020 


139319 


1871 


1674 


12 


147 


27 


31 


- 1021 


68074 


1024 


903 


7 


78 


17 


20 




192389 


()088 


5271 


i:! 


532 


109 


130 



TABLE I. Event numbers for the effect and estimated background at various energies. 
A^Stt - number of selected e+e" — + 7r"'"7r~7r° events after background subtraction. Ni,kg - 
the total estimated background. Contributions of separate background sources are also 
shown. The last line contains the total events numbers for all energy points. 



i 


rrii 
(MeV) 


"i,i-4, 


"i,i-3 


^(0) 








u(0) 
"i,i+2 


1,(0) 




hI°> ■ 10* 


1 


290 


- 










- 


0.001 


0.002 


- 


- 


1.055 


0.10 ±0.04 


2 


310 


- 










- 


0.015 


0.009 


0.001 


- 


1.004 


3.0 ±0.3 


3 


330 


- 










0.009 


0.210 


0.035 


0.001 


- 


0.018 


49.6 ± 1.1 


4 


350 




0, 


.001 


0, 


.002 


0.026 


0.295 


0.045 


0.002 


0.001 


0.015 


103.6 ±1.5 


5 


370 




0, 


.001 


0, 


.001 


0.028 


0.302 


0.051 


0.003 


0.001 


0.014 


154.0 ± 1.9 


6 


390 


0.001 





.001 





.003 


0.034 


0.317 


0.058 


0.003 


0.001 


0.013 


212.5 ±2.3 


7 


410 


0.001 


0, 


.001 





.002 


0.036 


0.316 


0.061 


0.003 


0.001 


0.012 


267.9 ± 2.5 


8 


430 











.001 


0.040 


0.313 


0.065 


0.004 


0.001 


0.012 


322.9 ± 2.8 


9 


450 











.002 


0.047 


0.304 


0.070 


0.004 


0.001 


0.012 


378.2 ±3.0 


10 


470 











.002 


0.050 


0.292 


0.071 


0.005 


0.001 


0.012 


434.7 ± 3.2 


11 


490 




0, 


.001 


0, 


.003 


0.058 


0.291 


0.075 


0.005 


0.001 


0.011 


485.9 ± 3.4 


12 


510 







.001 





.003 


0.061 


0.273 


0.073 


0.006 


0.001 


0.011 


539.0 ± 3.6 


13 


530 











.004 


0.065 


0.271 


0.076 


0.006 


0.001 


0.011 


578.3 ±3.8 


14 


550 




0, 


.001 


0, 


.005 


0.067 


0.258 


0.073 


0.006 


0.001 


0.011 


614.6 ± 3.9 


15 


570 




0, 


.001 


0, 


.005 


0.070 


0.249 


0.073 


0.006 


0.001 


0.011 


645.6 ± 4.0 


16 


590 




0, 


.001 


0, 


.006 


0.072 


0.239 


0.071 


0.006 


0.001 


0.011 


663.1 ±4.1 


17 


610 







.001 





.006 


0.071 


0.228 


0.069 


0.006 


0.001 


0.011 


668.3 ± 4.1 


18 


630 




0, 


.001 


0, 


.006 


0.070 


0.222 


0.067 


0.006 


0.001 


0.011 


650.9 ±4.0 


19 


650 







.001 


0, 


.006 


0.067 


0.210 


0.064 


0.005 




0.011 


628.8 ± 4.0 


20 


670 







.001 





.006 


0.065 


0.204 


0.059 


0.004 




0.011 


595.8 ± 4.0 


21 


690 







.001 





.006 


0.063 


0.195 


0.055 


0.003 




0.011 


544.5 ±3.8 


22 


710 




0, 


.001 





.005 


0.060 


0.187 


0.053 


0.003 




0.012 


472.5 ±3.7 


23 


730 











.005 


0.055 


0.179 


0.052 


0.002 




0.012 


380.7 ±3.5 


24 


750 








0, 


.003 


0.050 


0.169 


0.049 


0.002 




0.013 


287.2 ± 2.9 


25 


770 











.003 


0.045 


0.165 


0.042 


0.001 




0.014 


181.0 ±2.2 


26 


790 











.002 


0.035 


0.154 


0.034 






0.017 


90.1 ± 1.6 


27 


810 











.001 


0.026 


0.142 


0.027 






0.024 


32.3 ±0.9 


28 


830 












0.021 


0.113 


0.034 






0.039 


11.5 ±0.6 


29 


850 












0.010 


0.075 








1.004 


3.4 ± 0.4 


30 


870 












0.001 


0.045 








1.023 


0.10 ±0.06 



TABLE 11. The normalized ■w'^t: invariant mass spectrum, i/j - the relative fraction 
of events with 7r~^n~ invariant mass rrii ± 10 MeV, efficiency matrix b\°^ = a\°^ ■ e^p and 
relative uncertainty SP^^^ of the theoretical 7r"'"7r~ mass distribution are also given. 
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i 


rm 
(MeV) 




"i,i-3 


"i,i-2 




b 


if 












iff • 10^ 


1 


290 


— 




— 




— 


— 





030 


0.006 




— 




— 


_ 


1.004 


2.1 ±0.3 


2 


310 


— 




— 




— 


0.004 





081 


0.019 




- 




— 


_ 


1.001 


12.2 ±0.7 


3 


330 


— 




— 




- 


0.015 





119 


0.033 





001 




— 


_ 


0.021 


33.8 ± 1.1 


4 


350 


— 





001 





001 


0.020 





143 


0.046 





002 




— 


_ 


0.016 


62.0 ± 1.3 


5 


370 


— 




— 





001 


0.033 





162 


0.059 





004 




— 


_ 


0.013 


102.0 ± 1.7 


6 


390 


0.001 




— 





002 


0.039 





171 


0.064 





005 




— 


_ 


0.011 


144.9 ± 1.9 


7 


410 


_ 




— 





002 


0.048 





177 


0.071 





006 




- 


_ 


0.010 


194.4 ±2.3 


8 


430 


_ 




— 





002 


0.056 





183 


0.077 





008 





001 


_ 


0.009 


244.5 ± 2.6 


9 


450 


_ 




— 





004 


0.062 





185 


0.079 





010 





001 


_ 


0.008 


298.5 ± 2.8 


10 


470 


_ 




— 





005 


0.067 





190 


0.085 





Oil 





001 


_ 


0.007 


355.6 ± 3.0 


11 


490 


_ 





001 





007 


0.068 





185 


0.087 





013 





001 


0.001 


0.007 


417.9 ± 3.3 


12 


510 


_ 




- 





008 


0.074 





189 


0.088 





014 





002 


0.001 


0.006 


474.7 ±3.5 


13 


530 


_ 




— 





010 


0.075 





186 


0.091 





017 





003 


0.001 


0.006 


526.1 ±3.7 


14 


550 


_ 





001 





010 


0.076 





183 


0.090 





017 





003 


0.001 


0.006 


578.1 ± 3.9 


15 


570 


_ 





001 





013 


0.078 





177 


0.093 





017 





003 


0.001 


0.006 


620.4 ± 4.0 


16 


590 


_ 





001 





013 


0.081 





174 


0.094 





019 





004 


0.001 


0.006 


645.8 ± 4.1 


17 


610 


_ 





002 





015 


0.079 





170 


0.094 





019 





004 


0.001 


0.006 


667.3 ±4.1 


18 


630 


_ 





002 





015 


0.079 





168 


0.091 





020 





005 


0.001 


0.006 


676.0 ±4.1 


19 


650 


0.001 





003 





016 


0.080 





166 


0.090 





020 





005 


0.001 


0.006 


673.4 ±4.1 


20 


670 


0.001 





003 





016 


0.078 





163 


0.089 





021 





004 


0.001 


0.006 


660.6 ±4.1 


21 


690 


0.001 





003 





017 


0.078 





163 


0.090 





023 





005 


0.001 


0.006 


631.5 ±4.0 


22 


710 


0.001 





003 





016 


0.078 





162 


0.093 





021 





005 


0.001 


0.006 


570.7 ±3.8 


23 


730 


0.001 





004 





017 


0.078 





163 


0.090 





022 





005 


0.001 


0.007 


498.4 ± 3.5 


24 


750 


0.001 





004 





017 


0.077 





166 


0.094 





023 





004 


0.001 


0.008 


381.2 ± 3.0 


25 


770 


0.001 





003 





017 


0.077 





160 


0.091 





020 





002 




0.010 


259.5 ± 2.5 


26 


790 


0.001 





004 





017 


0.076 





171 


0.096 





018 





002 




0.012 


147.6 ± 1.9 


27 


810 


0.001 





003 





016 


0.073 





167 


0.087 





014 








0.018 


72.6 ± 1.3 


28 


830 


0.001 





003 





013 


0.070 





170 


0.054 












0.029 


30.8 ± 0.9 


29 


850 


0.001 





002 





013 


0.060 





171 


0.029 












1.001 


11.7 ±0.5 


30 


870 


0.001 





001 





012 


0.045 





057 














1.005 


5.6 ±0.4 



TABLE in. The normalized tt^tt" invariant mass spectrum, -ff}*' - the relative fraction 
of events with tt^tt" invariant mass rrii ± 10 MeV, efflciency matrix fo'^' = a^^' ■ e^*' and 
relative uncertainty (J-P/*' of the theoretical tt^tt*^ mass distribution are also given. 





1 


2 


3 


4 


m^o - 770 (MeV) 


6.1±1.0 


6.1±1.0 


4.9±1.4 


7.8±1.6 


Tpo (McV) 


151.1±2.6 


150.7±3.5 


156.0±4.0 


151.2±2.6 


fWp± — rripO (MeV) 
r^±(MeV) 


-1.7±1.2 


-1.8±1.2 


-0.6±1.6 


-1.9±1.2 


149.9±2.3 


149.5±3.4 


149.4±2.3 


150.2±2.3 


Re(o3^) • 10'^(MeV)-2 





0.04±0.25 
















0.3±0.1 













-90.0±18.0 





Si 











0.3±0.3 


xl 


27.1 


26.9 


19.0 


25.0 


xl 


28.9 


29.1 


29.5 


29.6 


xVN,f 


56.0/54 


56.0/53 


48.5/52 


54.6/53 



TABLE IV. Fit results for the tt^tt and tt^tt" mass spectra. 
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N 


1 


2 


3 


4 


rripo - 770 (MeV) 


5.8±0.9 


5.8±1.0 


4.3±1.3 


7.3±1.3 


rpo(MeV) 


149.9±2.2 


149.5±3.1 


149.3±2.3 


149.9±2.2 


rripi — rUpO (MeV) 
Re(a3^) • 10'5(MeV)-2 


-1.3±1.1 


-1.3±1.1 


0.3±1.5 


-1.4± 1.1 





0.04±0.24 
















0.2±0.1 













-125.U±2o.U 





Si 





U 


U 


n o_l_n o 




26.9 


26.8 


24.0 


24.9 


xl 


29.9 


oU.U 


oU.U 


oU. / 


xVNjf 


56.8/55 


OD.O/ 


0^.0/ Do 


oO.Of 0^ 


TABLE V. Fit results for the tt+tt" 


and tt^tt" 


mass spectra under 


assumption 












N 


1 


2 


3 


4 


rripo - 770 (MeV) 


5.0±0.6 


5.0±0.7 


4.5±0.7 


6.4±1.4 


rpo(MeV) 
Re(o3^) • 10'5(MeV)-2 


149.8±2.2 


149.7±3.1 


149.3±2.3 


149.9±2.2 





0.01±0.23 








-^(jjtt/ pir 








0.2±0.1 





^ujTv/ pn 








-125.0±28.0 





Sl 











0.3±0.3 


xl 


28.2 


28.0 


24.0 


26.2 


xl 


30.0 


30.0 


30.0 


30.8 


xVN,f 


58.2/56 


58.0/55 


58.0/54 


57.1/55 



TABLE VI. Fit results for the n n and n n mass spectra, under assumptions 
9p0nn = 9p±nn and rUpO = m^i. 
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mpo(MeV) 


m + fMeV") 


7?ip± — (JVIeV) 


SND 




775.8 ± 0.9 ± 2.0 


774.5 ± 0.7 ± 1.5 


— 1.3 ± 1.1 ± 2.0 


SND* r 




775 ± 6 ± 1 1 


775 ± 6 ± 1 1 




CMD-2il 


775.3 ± 0.6 ± 0.2 






ALEPH 


2(1 




776.4 ± 0.9 ± 1.5 


0.0 ± 1.0 


cLEok; 






774.9 ± 0.5 ± 0.9 




cbar| 


) 


762.3 ± 0.5 ± 1.2 


763.0 ± 0.3 ± 1.3 


-1.6 ± 0.6 ± 1.7 


PDG-lL 




776.0 ± 0.9 


776.0 ± 0.9 




PDG-21. 


2 


769.3 ±0.8 


769.3 ±0.8 


-0.4 ±0.8 






rpo(MeV) 


V(MeV) 


aa^ X 10"(MeV)"" 


SND 




151.1 ±2.6 ±3.0 


149.9 ±2.3 ±2.0 




SND**p 




149.8 ±2.2 ±2.0 


150.9 ±2.2 ±2.0 


0.01 ±0.23 ±0.25 


CMD-2I 




147.7 ± 1.3 ±0.4 




0.08 ±0.45 ±0.37 


ALEPJH 


20 




150.5 ± 1.6 ±6.3 




CLEOE 






149.0 ± 1.1 ±0.7 




cbarE. 




147.0 ±2.5 


149.5 ± 1.3 




PDG-lL 




150.5 ±2.7 


150.5 ±2.7 




PDG-2L 


2 


150.2 ±0.8 


150.2 ±0.8 












Sl 


SND 




0.20 ±0.10 ±0.05 


-125° ±28° 


0.3 ±0.3 ±0.3 



TABLE VII. Comparison of the results of various experiments. SND - this work; SND* 
- the p meson mass from our fit under m^o — rnp± assumption; SND** - the p° and 
widths from our fit under gpO„„ — gp±T^T^ assumption; PDG-1 - the world average values 
of e~ annihilation and r decay experiments only; PDG-2 - the world average values for 
all experiments. 
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7Z, 

y p,pl CO £^71 




FIG. 1. e"'"e~ ^ (j) ^ pn ^ -k'^-k~'k^ transition diagrams, a), b), c) - through various 
intermediate states, d) - the e^e~ ^ (f> ^ pn ^ tv'^tv'tt'^ transition with p and tt mesons 
interaction in the final state. 
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FIG. 4. The 6 distribution of neutral pions from the reaction e"'"e~ n^n~ 
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FIG. 5. The 6 distribution of charged pions from the reaction e"'"e~ — » 7r"'"7r~ 
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FIG. 9. "7r"'"7r~" mass distribution from the e^e~ — > KsKl background events. 
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FIG. 10. 7r"'"7r~ mass distribution for selected e^e~ — »■ 7r"'"7r~7r° events and background 
contribution 
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FIG. 11. 7r='=7r° mass distribution for selected e+e" — > 7r"'"7r~7r'' events and background 
contribution 
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FIG. 12. $(m, s) as the function of the dipion mass m aty'i = m,/, 
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FIG. 13. Detection efficiency for tv'^tt mass spectrum. 
^ 0.5 r 



E 



0.4 
0.3 
0.2 
0.1 














• 

• 






• 

• 

• 






i 

± 


• 

• 








• 

1 1 


1 1 1 


1 1 1 


Y' 



400 



600 800 

^0 (MeV) 



FIG. 14. Detection efficiency for tt n mass spectrum 
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FIG. 15. The correction coefficients p- ' and ' upto irrelevant normalization factor. 
Dots - correction coefficient to the detection efficiency for tv~^tv~ mass spectrum. 
Circles - /j'*^ correction coefficient to the detection efficiency for n^n° mass spectrum 




FIG. 16. The difference between the experimental and theoretical m^+^- spectra. 
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FIG. 17. The difference between the experimental and theoretical m^±^o spectra. 
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